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SUMMARY 


A method has been developed for the design and analysis of two- 
dimensional subsonic-flow passages with isentropic nonviscous flow. The 
method is based on the relation between the pressure change across a 
stream tube and the centrifugal force resulting from the curvature of the 
flow. Precomputed charts can be drawn which eliminate subsequent calcu- 
lation for a given upstream Mach number. The method is limited by the 
accuracy with which the radius of curvature of the streamlines can be 
determined . 

In an example the method was applied to the design of an expanding 
elbow at each of two Mach numbers . Two elbow contours and their surface 
pressure distributions were obtained. 


INTRODUCTION 


Much experimental work on internal-flow systems has consisted of 
tests of more or less arbitrarily laid out duct components to determine 
regions of separated flow, pressure losses, and flow distortion at the 
exit. This research could be conducted in a more effective manner if a 
method were available to design some of these duct components to avoid 
regions of local speedup and steep adverse pressure gradients. Further- 
more, if the flow in existing components could be analyzed, the analysis 
would furnish a basis for comparison of ideal and actual performance and 
could be used as a basis of elimination of unpromising components from 
an experimental program. 

Several methods for analyzing the flow in subsonic duct systems are 
available. Flugel's stream-filament method discussed in reference 1 was 
utilized in references 2 and 3 by making arbitrary assumptions about the 
variation of the radius of curvature . The method was used in reference k 
for the analysis of the flow in three-dimensional passages and in refer- 
ence 5 tor the design of three-dimensional passages. A wire-mesh plotting 
device for incompressible flow is described in reference 6, and in 


2 


NACA TN 4241 


reference 7 a device consisting of cams combined with a wire mesh is 
described which is applicable to compressible flow. 

In the present report a method is developed for the design and anal- 
ysis of two-dimensional subsonic-flow passages with isentropic nonviscous 
flow. The method parallels the three-dimensional treatment of refer- 
ences 4 and 5. Restriction of the treatment to two-dimensional-flow prob- 
lems permits the use of precomputed charts from which the values required 
can be read directly. The charts can be used for design or analysis of 
flow passages with the assumption of isentropic, subsonic, nonviscous 
flow. 


DEVELOPMENT OF METHOD 


The method developed in the present report is based on a relation 
between the pressure difference across a small stream tube and the cen- 
trifugal force per unit area caused by the flow curvature: 



In the following diagram is shown a small stream tube of width An 
and average radius of curvature r Q (all symbols are defined in appen- 
dix A): 



The ratio of the velocities at the intersection of a normal to the 
boundaries of the stream tube is derived in appendix B and given by 
equation (l): 

v An/?c. 

=t- = e 


(i) 
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Since no density values are involved, equation (l) is applicable to 
either incompressible flow or compressible flow. 

Continuity for two-dimensional flow requires that between any sta- 
tion and an upstream station where the flow is rectilinear 


pVAn = 


(2a) 


for compressible flow and 


VAn = 

for incompressible flow. 

If V is taken as V = i (V 0 + V), then 



(2b) 


(3) 


Velocity and density for isentropic compressible flow are related 
by equation (63) of reference 8: 



(4) 


By use of elementary isentropic relations and equations (2a), (3), 
and (4) it can be shown that 




V 




( 5 ) 


For consistency equation (l) is rewritten as 
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Equations ( 5 ) and (6) are obviously not easily solved for V/V u 
and An/Ar^ for given values of V 0 /V u , An u /? c , and M u . A chart can 
be made, however, to represent equations (5) and (6) for any given upstream 
Mach number by assigning a value to V Q /V U and using several values of 

^44 to obtain values of V/V , An/An u , and ^14 = ^14 4S_/4S_\ 

? c ^ ' U ' U r c ? c 

Thus lines for constant values of V Q /V U can be drawn on a chart with 

An/Ax^ for ordinate and An u /r Q for abscissa. By cross-plotting, lines 

of constant values of V/V u can be superposed on this chart. However, a 

simpler method is described in appendix C whereby a chart can be easily 
constructed for a given upstream Mach number. From this chart values of 
An / /An u and V /V u can be read directly for given values of V Q /V U 

and An u /r c . Figure 1 shows charts for incompressible flow. Figures 2 

and 3 show charts for compressible flow for upstream Mach numbers of 0.40 
and 0.85, respectively. If a value of An u can be decided upon for a 

particular case, the corresponding curves can be relabeled with values of 
An and l/r c . In figure 4, which is for an upstream value of An of 

2.00 inches, the values of An are given in inches and values of 1 /r c 

in (inches) "■4. The lines of this figure are the same as those of fig- 
ure 3 but have different labels . 

The accuracy of this method depends in part on the accuracy with 
which l/? c can be determined and decreases as l/r c becomes very large 

or very small. Discontinuities in r c for the starting boundary should 
be kept to a minimum magnitude. The correct sign for 1 jf c must be 
used. The sign of 1 /r c is the same as that of the second derivative 

of the curve and can usually be determined by inspection. If the curve 
is concave when viewed from the direction from the boundary to stream- 
line 1, the reciprocal of the radius of curvature is negative. Inadvert- 
ent use of the wrong sign will immediately be apparent because the direc- 
tion of pressure change indicated by the change in velocity across a 
stream tube will not be consistent with the direction of turning. It is 
considered that the accuracy decreases for successive streamlines because 
the measured values of l/r G depend more on the artistry of fairing in 
the streamlines through the An/Ar^ points. 
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APPLICATION OF MEOHOD 


Design of Duct Shape 

Use of the charts for design purposes is essentially as follows: 

A boundary or other streamline along which a chosen ^velocity variation 
is prescribed is laid out and perpendiculars are erebted at a suitable 
number of stations. If the boundary has been determined analytically, 
it may be possible to compute the radius of curvature at each station 
and estimate the average radius of curvature for the_ stream tube. Other- 
wise, some other means, such as the radometer shown 'in figure 5 and 
described in appendix D, is required to determine the radius of curvature. 
For each station, values of An/An u and V/V u are read at the proper 

values of V Q /V U and An u /r c - The first streamline can now be drawn in 

and the average radius of curvature determined by measurements. The new 
values of An u /r c based on the measured value of l/r c must be used to 
determine new values of An/Ar^ and V/V u . This iteration procedure must 
continue until the values of An/An^ are essentially unchanged by succes- 
sive values of An u /r c obtained from measurements. The number of repeti- 
tions required is reduced by the fact that changes in the values of An u /? c 
do not give large changes in the values of An/An u . For example, a 
50-percent increase in An u /? c (from 0.030 to 0.045) in figure 3(c) at 
V Q /V U of 0.99 gives a change in An/Ar^ of only 0.5 percent. 

The final values of V/V u are now used as values of V Q /V U for the 
next streamline. In constructing the next streamline the perpendicular 
along which An/An^ will be measured is laid off from the midpoint of 

the perpendicular used in the previous stream tube. The values of V/V u 

for the final streamline can by recourse to tables for compressible flow 
be converted to obtain pres sure -distribution curves. 


Analysis of Flow Quantities 

The method of analysis is used when the inner and outer boundaries 
are prescribed. The procedure to be described is given as one possibility, 
but may, of course, be altered to suit the individual situation. 

After stations are selected along one of the boundaries, an approxi- 
mation to a line perpendicular to the streamlines must be drawn in for 
each station. Drawing in estimated streamlines might be of some assist- 
ance in making this approximation. Selection of the number of stream 
tubes to be used determines the upstream value of An. The total length 
of the line perpendicular to the streamlines is measured for each station. 
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At each station average values of An u /r c must be determined for 
each stream tube either analytically or by measuring the values of An u /r c 
at both boundaries and plotting a reasonable curve between the two values. 

With the streamlines numbered from 0 for the starting boundary up to 
the other boundary, which will have a number equal to the number of stream 
tubes, a tabulation can be made for each station as shown in the fol- 
lowing table: 


Stream 

tube 

number 

^u 

? c 

v 0 

An 

An u 

V o 

v u 

An 

Anu 

1 

2 

3 






Summation of values 

An u 


In the table the stream tubes are identified with the streamline whose 
location is being determined. The values of An u /r c already determined 

are written in for each stream tube. 

A value of V Q /V U is needed for the first boundary at each station. 

It can be obtained from the precomputed charts at the proper value of 
An u /r c and an estimated value of An/An u . This value of V Q /V U is 

written in the first column headed V Q /V U for stream tube 1. The pre- 
computed chart is used to determine values of An/An u and V/V u for 
stream tube 1; however, the value of V /V u for stream tube 1 is written 
in as the value of V Q /V U for stream tube 2. This procedure is con- 
tinued until there is a value of An/An u for each stream tube and one 

more velocity ratio than there are stream tubes. The sum of the computed 
values of An/An^ is compared with the sum of the values of An/An u 

measured on the layout. If these sums are not equal, a new value of 
V 0 /V u for stream tube 1 is inserted in the second column headed V 0 /V u 

and new values are picked from the precomputed chart. The new initial 
value of V Q /V U is obtained by multiplying the value previously assumed 

by the ratio of the sum of the computed values of An/An^ to the sum of 
the measured values of An/Ar^. This routine may have to be repeated if 
the sums of An/An u are still not sufficiently close. 
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The computed values of An/Ar^ are now used to plot streamlines. 

For each station, perpendiculars to the streamlines are constructed in 
the manner described for the design method. Values of An/r c determined 

from measurements are inserted in the table and a value of the summation 
of An/Ar^ is determined from measured values of An/Ar^. A value of 
V 0 /V u f° r the first boundary is obtained from the precomputed charts 
and values of An/Ar^ and V/V u are successively obtained for each 
stream tube. This procedure is repeated with a new value of V 0 /V u for 
the start if the sum of the computed values of An /An^ does not equal 

the sum of the measured values . If necessary, new streamlines are drawn 
and the process repeated until the streamlines plotted from computed 
values of An /Ar^ are essentially unchanged because of the use of meas- 
ured values of An u /r c . The required number of repetitions is reduced 
by the fact that, as has been pointed out, changes in An u /r c do not 
cause proportionate changes in An/Ar^. 


Examples 

The design method was used to lay out an expanding 90 ° elbow 
(fig. 6) for an upstream Mach number of O.85. The area ratio was l- 37 > 
and the space used up was to be confined to moderate bulges outside that 
required by a nonexpanding elbow of mean radius three times the inlet 
width. It was thought that use of an entering Mach number higher than 
O.85 might result in some supersonic flow in the elbow. 

A modified sine curve was used for the inner boundary in order to 
have zero values of l/r c at the beginning and the end of the turn. It 
was assumed advantageous to spread the adverse pressure gradient evenly 
over the entire length of the inner surface. 

For the layout an upstream width of 10 inches was divided into 
five stream tubes having equal flow. 

It was found that some drop in pressure on the inner wall in the 
region of the start of the turn was required in order to avoid a necking-in 
of the outer wall. This necessity for a drop in pressure can be explained 
as follows: Consider first a station a short distance past the start of 

the turning of the inner surface. The assumed distributed pressure rise 
gives a small increase in pressure over the upstream value. The turning 
requires that the next streamline out have a slightly greater increase in 
pressure and that, consequently, the stream-tube width be increased over 
the upstream value. By this same reasoning each successive stream- tube 
width at this station will be slightly greater than that of the preceding 
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one. These increases in stream- tube width move the outer boundary at 
this station outside of a line which is an extension of the outer straight 
inlet boundary. This outward movement of the boundary means that at some 
previous station the outer boundary had to curve outward to reach this 
point. Now consider this previous station. At the inner boundary there 
is no curvature and the pressure has not changed from the upstream value. 
For the next streamlines the pressure and stream-tube width are unchanged 
up to the region where the outer streamlines are curved outward; there- 
fore, each of the last streamlines must have a lower pressure than the 
one before it and the stream- tube widths are successively decreasing. 

The result is that the outer boundary is pulled inside the line of the 
straight inlet. It was assumed that a decrease in area ahead of a dif- 
fusing bend was not favorable . The decrease in area, was eliminated by 
applying an inner boundary pressure lower than the upstream pressure to 
both of the stations previously discussed. At the station beyond the 
start of the bend this lower pressure resulted in moving all the stream- 
lines inward and reduced the amount that the outer boundary had been out- 
side the inlet line. The outward curvature of the outer boundary at the 
preceding station was therefore reduced. The application of the lower 
inner boundary pressure and the resultant moving inward of the stream- 
lines gave a curvature in the direction of the bend to the first stream- 
lines of the station previous to the start of the inner-wall bend. The 
result at this station was that the stream tubes near the inner walls 
were thinner than the upstream values. This result permitted the stream 
tubes in the center to be wider and the pressures to be higher and con- 
sistent with curvature away from the center of the passage toward each 
boundary. The prescribed inner -wall pressure was modified several times 
and the design was restarted until the result was an outer boundary which 
did not reduce the cross-section width before the start of the bend. 

The resultant elbow and streamlines are shown as the solid lines in 
figure 6 . The variation of velocity, stream-tube width, and radius of 
curvature along the streamlines is indicated in figure 7 . The corre- 
sponding pressure distributions over the inner and outer surfaces are 
shown as the solid lines in figure 8 . 

For comparison with this elbow designed for a Mach number of O .85 
the design of an elbow for an upstream Mach number of 0.40 was started. 

The Mach number of 0.40 was chosen as being in the region of effects 
midway between the example Mach number of O .85 and incompressible flow. 

The same inner-wall shape and area ratio were used. It was found that 
a larger drop in pressure at the start had to be tolerated. This drop 
in pressure is due in part to the fact that a greater area change is 
needed at the lower Mach number for a corresponding change in pressure. 

Toward the outlet of the elbow it became very difficult to modify 
the inner-wall velocity distribution in such a way as to get the outer 
boundary to come out smooth and fair with the outlet line and the first 
half of the outer boundary. The first half of the outer boundary was. 
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- therefore, faired into a curve arbitrarily taken from the last part of 

\ the outer boundary of the elbow designed for a Mach number of O.85, and 

> the analysis method was applied. 

\ The use of the analysis method permitted the intermediate stream- 

lines and the pressure distributions over the inner and outer surface 
to be obtained. The boundaries and streamlines are shown as the dashed 
lines in figure 6. The surface pressure distributions are shown as the 
dashed lines in figure 8. 

In the examples a 10-inch inlet and a value of An^ of 2.00 inches 

were used. It was found convenient to use charts similar to figure 4, 
which was made for a value of An^ of 2.00 inches, and to use a radometer 

which gave direct readings of l/r c in (inches) - ^. It was found helpful 
to make plots of l/r c and An against linear distance along the bound- 
ary for the design method. For the analytical method plots were made of 
l/r c and An against n. 

A rough check was made on the effects of the number of stream tubes 
» into which the passage is divided. Station 11 of the example for an 

upstream Mach number of 0.40 was investigated by the analytical method 
with 10 stream tubes instead of 5 . The solid line in figure 9 gives the 
pressure variation across the passage with the use of 5 stream tubes. 

The dashed line in figure 9 gives the pressure variation for 10 stream 
tubes. The greatest difference was O.OlOq at the inner boundary. Con- 
sistency would indicate an approximately equal change at the adjacent sta- 
tions so that the effect on the pressure gradient would be nil. The 
largest shift in streamline location was 0.02 inch. In this case the 
extra work required to use 10 stream tubes instead of 5 would apparently 
not be justified. 

CONCLUDING REMARKS 


In the present report a method has been developed for the design 
and analysis of two-dimensional subsonic-flow passages with isentropic 
nonviscous flow. The method is based on the relation between the pres- 
sure change across a stream tube and the centrifugal force resulting from 
the curvature of the flow. Precomputed charts can be drawn which eliminate 
subsequent calculation for a given upstream Mach number. The method is 
limited by the accuracy with which the radius of curvature of the stream- 
lines can be determined. 
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In an example the method was applied to the design of an expanding 
elbow at each of two Mach numbers . Two elbow contours and their surface 
pressure distributions were obtained. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 15, 1958. 
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APPENDIX A 
SYMBOLS 


a speed of sound, ft/sec 

g acceleration due to gravity, ft/sec^ 

M Mach number 

n distance along normal to streamline 

Ztfi distance between adjacent streamlines along normal 

p static pressure, lb/sq ft 

q dynamic pressure, ^ V^, lb/sq ft 


r c radius of curvature 

V velocity, ft/sec 

V velocity at first streamline bounding a stream tube, ft/sec 
o 

p / I}- 

p mass density, lb-sec^/ft 

y ratio of specific heats 

Subscripts : 


t stagnation conditions 

u conditions upstream where flow is not turning 

0,1, . . . points at succeeding streamlines along normal 

A bar over a symbol denotes an average value representative of the 
stream tube . 
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APPENDIX B 

DERIVATION OF VELOCITY EQUATION 


The velocity equation used in the method of the present report is 
derived from a relation between the pressure difference across a small 
stream tube and the centrifugal force per unit area caused by the flow 
curvature : 


or 


dp = -p 


dn — 
r c 


P 



From Euler's equation for steady one-dimensional flow. 


ftp 

P 


= -V dV 


Combining equations (Bl) and (B2) gives 


ftp 


P 


-V 



(Bl) 


(B2) 


Dividing by V^ 


yields 


d V = 
V 




ioge 
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l 


dn 


If An, 
value of r c 


the distance between the streamlines, 
may be used and, 


is small, an average 


V 


An/r„ 
V„e ' c 


(B3) 


Equation (B3) could have been obtained from equation (3) of refer- 
ence 5 by eliminating the last term of equation (3) which is zero in the 
case of two-dimensional flow. 
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APPENDIX C 

CONSTRUCTION OF PRECOMPUTED CHART FOR 
SELECTED UPSTREAM MACH NUMBER 


For a given upstream Mach number a chart can be constructed from 
which values of An/An^ and V/V u can be read directly for given values 

of V 0 /V u and An u /r c . In a plot of V/V u against V 0 /V u , a line 
through the origin with a slope of 1 represents the case of no curvature, 
— = 0. For incompressible flow (fig. 1), equation (2b) can be used to 

r c 

spot values of AnyAn u on this line. For compressible flow at a given 

upstream Mach number (figs. 2 and 3) tables for one-dimensional flow 
(ref. 8) can be used to determine values of An/Ar^. Since An/Ar^ 
determines V/V u for a fixed upstream Mach number, the averaging equation 



shows that lines drawn through the points on the line for — = 0 with a 

f c 

slope of -1 each represent conditions for a constant value of An/Ar^. 


Points spotted on several of the lines of constant values of An /An u 
for a chosen value of An^/r c permit drawing in a line for the constant 
value of An u /r c . These values are determined as follows: For a constant 

value of An/A it can be shown that the average velocity V/V u is 

. Use of this fact in com- 

= 0 

bination with equations (3) and (6) for a constant value of An/An^ 
gives 


equal to the no-curvature velocity I — 


u/ 


Yo 

V u 
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The change in V Q /V U due to curvature when An/An u is held constant 
is then 







0 



SI.. 

? c ^ r c 


An u An 

e - 1 

^u Ar^ 


(Cl) 


Three or four values of An suffice to establish a line for a given 
value of An u /r c since it is nearly a straight line. The complication 
of using equation (C 1) can be greatly reduced by use of the fact that the 


term 


An 

e r c ^ 

A% An 
r„ ^ 


is closely approximated by i ^ 

2 r c Anu 


e ** +1 

small ; therefore , 


^^7 lf An/r c 


is 



1 AHu An ( Y o 

2 r c 



(C2) 


Use of equation (C2) permits determination of the change in V Q /V U due 
to curvature with sufficient accuracy by simple multiplication. 
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For incompressible flow 



= 1 


and 

V o ^ 1 
V u 2 r c 

Lines of constant values of An u /r c for incompressible flow should then 

be a constant distance from the no-curvature line regardless of the value 
of An /Ar^ . 

Since this treatment is for subsonic flow neither V Q nor V may 

exceed the velocity of sound. The dashed limit lines in figures 3(c) 
and 4(c) for a value of M u of 0.85 represent this limitation. 

In the derivation (appendix B), n and r c would be expressed in 

feet to be consistent with the basic equations. In the derived equa- 
tion (l) and in the continuity equation (2) it will be noted that An 
appears only as a ratio to another An value or as a ratio to r c . In 

this case it is only necessary that An and r c be expressed in the 
same units of length. 
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APPENDIX D 

DETERMIMTION OF AN AVERAGE RADIUS OF CURVATURE 

A radometer is one means of determining a value of l/r c of a 

curve for which this value is not otherwise known. The radometer of 
figure 5 I s similar to the one described in reference 9* It gives, how- 
ever, a direct reading of l/r c . In figure 5(b), AC is a curve for 

which the value of l/r c is 1° be determined at point B. Points D 

a nd E equidistant from B are assumed to be close enough to B for this 
curve to exclude excessive change in l/r c between them. 

In the nomenclature of the figure, and Since 

bd can be a constant for a given radometer, the arc with radius b and 
center at B can be labeled to read values of l/r c . 

In use, the pivot point B is placed on the curve at the point at 
which l/r c is to be estimated, the point D of the stationary part is 
put on the curve, and the movable part is moved until the point E also 
is on the curve. The value of l/r c can now be read. 

For the instrument used for the example of this report, b was 

10 inches and d was 5 inches; l/r c is expressed in (inches) -1 for 
this radometer. 
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(a) 0.50 < — < 1.00. 

V u 


Figure 1.- Design chart. Incompressible flow 
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(b) o.8o < Ja < 1.30. 

v u 


Figure 1.- Continued. 
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Figure 1.- Concluded. 
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Figure 2.- Continued. 
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(c) 1.10 < ^ < 1.60. 

V u 


Figure 2.- Concluded. 
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Figure 3.- Design chart. Upstream Mach number, O.85. 



> > 


26 


NACA TN lj-241 



(b) 0.50 < ^ < 1.00. 

v u 


Figure 3-- Continued. 


< c 


NACA TN 4241 


27 



V 

(c) 0.70 < < 1.20. 

v u 


Figure 3*- Concluded. 
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Figure 4.- Specific design chart for upstream width An^ of 2.00 inches 

Upstream Mach number, 0.85* 
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Figure 4.- Continued 
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(c) 0.70 < < 1.20. 

v u 


Figure 4.- Concluded. 
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(a) Photograph of instrument. L-57-2564 


Figure 5«- Radometer. 
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(b) Schematic diagram. 


Figure 5-- Concluded. 
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Figure 6.- Elbow contours and streamlines resulting from an arbitrarily 
assigned inner boundary and an outlet area 1.37 times inlet area. 
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Station number 


(a) Variation of 

Figure 7 .- Variation of velocity and geometric quantities in a 90° diffusing elbow with an 

upstream Mach number of O. 85 . 
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(b) Variation of An. 


Figure 7 .- Continued. 
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(c) Variation of — . 

r c 


Figure 7*- Concluded. 
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Figure 8.- 


Comparison of pressure distribution for an upstream Mach number of 0.40 with that 
for an upstream Mach number of 0.85 . 
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Figure 9*- Comparison of pressure distribution determined for 10 stream tubes with that deter- 
mined for 5 stream tubes. Upstream Mach number, 0.-40. 
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